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Abstract: Based on the hard-thermal-loop resummed improved ladder Dyson- Schwinger 
'[^" ' equation for the fermion mass function, we propose a procedure how we can get the gauge 
O I invariant solution in the sense it satisfies the Ward-Takahashi identity. Results of the numerical 
^ I analysis are shown and properties of the "gauge-invariant" solutions are discussed. 

1. Introduction and Summary 

^ , 

The Dyson- Schwinger equation (DSE) has proven itself a powerful tool to investigate with 
' r-| '! the analytic procedure the nonperturbative structure of field theories, such as the chiral phase 
Q-i| transition of gauge theories. Actually analyses based on the DSE have successfully clarified 
Qh| the phase structure of vacuum QEC/QCD [1,2,3]. Here we must take note of the fact that 
^ ■ these DSE analyses in vacuum gauge theories were carried out in the Landau gauge with the 
ladder interaction kernel. Same analyses have been performed at finite temperature(T) and/or 
density(/i) also in the Landau gauge with the ladder interaction kernel [4,5,6], and with the 
hard-thermal-loop (HTL) resummed improved ladder interaction kernel [7]. 
Cs) ! The reason why the ladder approximation is used is that the full DSEs are coupled inte- 

^ I gral equations for several unknown functions, thus are hard to be solved without introducing 
I appropriate approximations. We usually adopt the step-by-step approach to this problem, 
O ' firstly approximate the integration kernel by the tree, or, ladder interaction kernel, next use 



> 



o 



the improved ladder one, etc. The possibility of such a systematic improvement through the 

well-established analytic procedure is one of the important characteristics of the DSE. In fact, 

. results of Ref.[7] show that at finite temperature/density it is important to correctly take the 

! dominant thermal effect into the interaction kernel in terms of the HTL resummation. 
H , 

. 5t , In case of the vacuum QED in the Landau gauge DSE with the ladder kernel for the fermion 

mass function, the fermion wave function renormalization constant is guaranteed to be unity 
[1], namely the Ward-Takahashi identity is automatically satisfied. Thus irrespective of the 
problem of the ladder approximation, the results obtained would be gauge invariant. 

At finite T and/or /i, however, there is no such guarantee. In fact, even in the Landau gauge 
the fermion wave function renormalization constant largely deviates from unity [7,8], being not 
even real. At finite T and/or /i the results obtained from the ladder DSE explicitly violate 
the Ward-Takahashi identity, thus depend on the choice of gauge. All the preceding analyses 
[4-7] suffer from the gauge-dependence problem coming from the ladder approximation of the 
interaction kernel, their physical meaning being obscure. 

In this paper, we present, in the analysis of the HTL resummed improved ladder DS equation 
for the fermion mass function in thermal QED/QCD, the procedure to get the gauge invariant 
solution in the sense it satisfies the Ward-Takahashi identity. We firstly show that the solutions 
of the HTL resummed improved ladder DS equation in thermal QED/QCD suffer from the 
problem of gauge-parameter dependence, then solve numerically the DSE constrained to satisfy 
the Ward-Takahashi identity and investigate the properties of the "gauge invariant" solution. 



Part of the preliminary result of the analysis was reported in Ref. [9] , showing the effectiveness 
of the procedure. 

Results of the present analysis are summarized as follows: 

(1) The solution of the HTL resummed improved ladder DS equation depends strongly on the 
choice of the gauge parameter within the momentum-independent gauge. This type of solution 
always shows the explicit contradiction with the Ward-Takahashi identity. 

(2) We can determine numerically the solution that satisfies the Ward-Takahashi identity, 
namely the solution in which the fermion wave function renormalization constant is almost 
equal to unity. To get such a solution it is essential that we work in the nonlinear gauge where 
the gauge parameter ^ depends on the momentum of the gauge boson. 

(3) The chiral phase transition in the massless thermal QED/QCD is confirmed to occur through 
the second order transition; a dynamical fermion mass is generated at the critical temperature 
or at the critical coupling constant without discontinuity. 

(4) The effect of thermal fluctuation on the chiral symmetry breaking and/or restoration is 
smaller than that expected in the previous analysis in the Landau gauge [7]. 

2. DS equation for the fermion self-energy function 

2-1. DS equation in the HTL resummed improved ladder approximation 

The fermion self-energy function appearing in the fermion propagator Sr 

Sn{P)^[r+ie^'-^RiPr' (1) 

can be decomposed at finite temperature and/or density as 

E«(P) = (1 - A{P))pa' - B{P)j' + C{P) (2) 

with A[P), B[P) and C(P) being the three scalar invariants to be determined. In the present 
analysis, we use the HTL resummed form *G^i, for the gauge boson propagator, 

1 ^ 1 V ^ V 

where *IIl/t is the HTL resummed longitudinal/transverse photon self-energy function [10]. 
A'^", B^^ and D^^ are the projection tensors [11], 

^M^^ ^ gi^^ _ Bi''' - D"'' , (4) 
Bi^i^ = -Ki'K''/K'^, (5) 
jji^i^ = K^K'^jK'^, (6) 

where K = (k, fcok), k = \^ and k = k/A; denotes the unit three vector along k. 

The parameter ^ appearing in the term proportional to the projection tensor D^i, represents 
the gauge-fixing parameter = in the Landau gauge). This gauge term plays an important 
role in the present analysis. 



The vertex function is approximated by the tree (point) vertex. With the instantaneous 
exchange approximation for the longitudinal gauge boson propagator, we get the DSEs for the 
three invariant functions A(P), B(P) and C(P) 
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\ko + B{K)+ief - A{KYk'' -C{Kf^\ ' 

Above DSEs may have several solutions, and we choose the "true" solution by evaluating 
the effective potential V^[S'ij] for the fermion propagator function Sr, then finding the lowest 
energy solution. 

V[Sr] = zIV[f5fl]+zTrlntei 



e 



d^K r d^P 1 
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+j^SciK)j.Sn{P)D^^^{P -K)+ j,Sr{K)j,Sc{P)D^/{P - K)] , (10) 



2-2. Procedure to get the "gauge-invariant" solution 

The function A{P) above is nothing but the inverse of the fermion wave function renormal- 
ization constant Z2, thus must be unity in order to satisfy the Ward-Takahashi identity in the 
ladder DSE analysis, where the vertex function receives no renormalization effect, Zi = 1. 

We must solve the above DSEs and get the solution satisfying the Ward-Takahashi identity 
Z2 — Z\{— 1), where Z2 — A{P)~^ . The procedure to get the "gauge invariant" solution is as 
follows; 

(1) Assume the nonlinear gauge such that the gauge parameter ^ being a function of the 
momentum K = {ko, k) carried by the gauge boson. We parametrize ^ as 

C{ko,k)^Y.^mnHm{ko)Ln{k), k ^ V¥ , (11) 

where ,^„j„ are unknown parameters to be determined. Hm and L„ can in general be any ortho- 
normal functions, and we here take the Hermite functions for and the Laguerre functions 
for L„. 

(2) When solving the above DSEs iteratively, impose the condition A{P) = 1 by constraint for 
the input-functions at each step of the iteration. 

(3) Determine so as to minimize \ A{P) — Ip for the output-functions and find the solutions 
for B{P) and C{P). 

3. "Gauge invariant" solution consistent with the Ward-Takahashi 
identity 

Here we present the results obtained by the momentum-dependent guge parameter ^. Num- 
ber of parameters ^mn to minimize \A{P) — Ip is 2 x 3 x 2 = 12 (i.e., m = ~ 2 and n — 0,1) 
in the case of complex ^, and 4 x 3 = 12 (i.e., m = --^ 3 and n = ~ 2) in the case of real 
All the quantities with the mass dimension are evaluated in the unit of A, the cut-off parameter 
introduced as usual to regularize the DSEs. 

Before presenting the "gauge invariant" solution, we show in Fig.l the result of the critical 
temperature analysis for several values of constant ^ to get a rough image for the size of gauge 
dependence. As can be seen clearly the critical temperature strongly depends on the gauge, 
but the order of the phase transition docs not. 

Now we present the solution consistent with the Ward-Takahashi identity, i.e., the "gauge 
invariant" solution. Analysis is now in progress, and the results shown below are, at present, 
still preliminary. 

Firstly in Fig.2 we show Re[A{P)]. For comparison, results in the constant ^ analyses are 
also shown in the same figure. 

Next let us study the property of the phase transition. Fig.3 shows the real part of the 
fermion mass Re[M{P)], M{P) = C{P)/A{P) (= C(P), because A{P) = 1), obtained from 
the "gauge invariant" solution, as a function of the temperature T. The mass is evaluated at 
Po = 0, p = 0.1, to be consistent with the standard prescription to define the mass in the static 
limit, Po = 0, p — > 0. 
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Figure 1: Gauge-par ameter-depnendence of the fermion mass Re[M] at the couphng constant a = 4.0 
evaluated at = 0, p = 0.1. 
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Figure 2: Comparison of the wave function renormalization constant Re [A] at the coupling constant 
a = 4.0 evaluated at ^>o = 0, p = 0.1. 
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Figure 3: Temperature-dependence of the fermion mass Re[M] for various values of tire coupling 
constant a = evaluated at po = 0, p = 0.1, see text. Open circle denotes real ^ data, while solid circle 
denotes complex ^ data. 

Analyses to determine the critical temperature Tc, the critical coupling etc, and two critical 
exponents u and r) are now in progress. 

As can be clearly seen in the above Figures 2 and 3, two solutions obtained in the different 
two prescriptions, complex v.s. real gauge parameters ^, show complete agreement. This fact 
indicates that the solution obtained in the present procedure does not depend on the choice of 
the gauge parameter, namely that the solution is "gauge invariant" The results shown in this 
paper are still preliminary, and we will soon report the results of full analysis. 

The phase boundary curve in the (T, Q;)-plane thus determined shows that the region of the 
symmetry broken phase shrinks to the low-temperature and the strong-coupling side compared 
with that of the Landau gauge. This fact means that the effect of thermal fluctuation on the 
chiral symmetry breaking/restoration is smaller than that expected in the previous analysis in 
the Landau gauge [7]. 

4. Discussion and comments 

Results presented in the present paper are still preliminary, because of the rough analysis 
of the data processing. We are now refining the data analysis and soon get the results of 
the thorough reanalysis. Though the main conclusion will not be altered, several important 
remarks should be added. 

(1) We performed the present analysis in two prescriptions for the nonlinear gauge parameter 
^, complex v.s. real. Allowing the gauge parameter ^ to be a complex value may correspond 
to studying the non-hermite dynamics, thus may cause some troubles. In this sense we are 
interested in the result obtained by restricting the gauge parameter to the real value. What 
we found is a remarkable result: In both cases results completely agree, thus getting a solution 
totally independent of the choice of gauges. This fact strongly indicates that we can get the 
gauge-invariant physical result by studying the DSE with the ladder interaction kernel through 



the present procedure. 

(2) In the present analysis, the consistency of the solution with the Ward-Takahashi identity is 
respected only by imposing the condition A{P) ^ 1. Needless to say, in solving the (improved) 
ladder Dyson- Schwinger equation, there are no solutions totally consistent with the Ward- 
Takahashi identity. Despite this fact, following point should be closely examined: At least 
around or in the static limit, po = 0, p ^ 0, where we calculated (defined) the mass, each 
invariant function A{P), B{P) or C{P) should not have big momentum dependence. This 
condition may be important in connection with the consistency of the obtained solution with 
the gauge invariance. Result of the present analysis shows that at least B{P) and C{P) satisfy 
this condition. 
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